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Abstract 

Qualitatively new transport mechanism is suggested for hopping of carriers 
according to which the variable-range hopping (VRH) arises from the reso- 
nant tunneling between transport states brought into resonance by Coulomb 
potentials produced by surrounding sites with fluctuating occupations. A 
semiquantitative description of the hopping transport is given based on the 
assumption that fluctuations of energies of hopping sites have spectral density 

V/- 



Recently, it was suggested that fluctuations of site occupations in ther- 
mal equilibrium can cause essential energy fluctuations for sites participat- 
ing in the hopping transport [1,2]. This idea was used to explain the low- 
frequency noise with spectral density 1/f reported for hopping conductiv- 
ity in p-Ge(Ga) and n-Ge(As) [1]. Due to the long-range character of the 
Coulomb potential, fluctuations in occupation numbers of some sites lead 
to essential fluctuations of energies of the surrounding sites. In Refs.[l,2], 
an attempt has been done to distinguish between the sites participating in 
the hopping transport (transport sites) and fluctuators, the latter being the 
sites which produce energy fluctuations on the transport sites, but do not 
directly participate in hopping processes. Comparison of experimental data 
with simple theoretical estimates [2] shows that the noise level is so large that 
each effective resistor on the transport cluster has relative fluctuations of the 
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resistivity of the order unity. In the present report we study the influence of 
such fluctuations on the variable-range hopping. 

According to the theoretical picture suggested in Refs.[l,2], the current 
noise is due to fluctuations of site energies affecting the hopping probabilities. 
The fluctuations of site energies are caused by electron transitions between 
sites which do not belong to the percolation cluster, but can influence the 
sites on the cluster by their Coulomb potentials. Thus this noise is a direct 
manifestation of the sequential Coulomb correlations in a hopping system 
suggested first by Pollak and Knotek [3,4]. The effect of these correlations 
on hopping transport has been discussed in numerous papers, although most 
studies in this field were based so far on computer simulations. Among 
very interesting results one can mention those of the recent study of Perez- 
Garrido et al. [5], who showed that the sequential correlations along with 
electron polarons have significant influence on the hopping conductivity in the 
Coulomb gap regime. Much less number of attempts have been made yet to 
develop analytical approaches to study the influence of sequential correlations 
on hopping conductivity. Ortuno and Pollak [6] were the first who considered 
analytically the influence of these correlations on the activation energy of 
hopping conductivity under the assumption of a sharp Coulomb gap. Using 
the mean-field theory they have shown that sequential correlations can reduce 
the width of the Coulomb gap Ac and hence diminish the activation energy 
of hopping conductivity. 

In the present report, we study analytically the effect of sequential correla- 
tions within a somewhat different approach. Following the ideas of Refs.[l,2], 
we presume, that sequential correlations cause fluctuations of site energies 
of the effective sites on the percolation cluster with spectral density 1//. In 
general case it is difficult, if possible at all, to distinguish between transport 
sites and fluctuators [7] and we are not able to perform this distinguishing in 
the whole frequency range. However, in the limit of the high-frequency noise 
it is reasonable to assume [1,2] that fluctuators are represented by compact 
pairs of sites each pair having one electron with correspondingly high transi- 
tion rate between the sites of the pair. Such pairs play the role of electronic 
two-level systems [2] and, as the well-known two level systems in structural 
glasses, lead to 1/ f spectrum of fluctuations. Thus in this frequency region 
the fluctuations of energies of hopping sites are caused by hopping ransitions 
of electrons within separated low-energy pairs surrounding the hopping sites 
similar to those considered by Efros et al. [8]. Below we will assume that 
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this 1/f fluctuations spectrum starting from the cut-off frequencies within 
the high-frequency region mentioned above is extended (with the same pa- 
rameters) to much lower frequencies where the fluctuations are dominated 
by the aggregates of defects more complex than the isolated pairs mentioned 
above. Having information neither about exact picture of these aggregates 
nor about the character of the corresponding relaxational modes, we still 
believe that they support the same 1/f fluctuation spectrum started at high 
frequencies where it is supported by isolated pairs. We would like to empha- 
sise that this assumption is decisive for the whole consideration suggested 
below. We cannot justify this assumption for the broad frequency range and 
we take it just as an ansatz. We also assume the fluctuations of site energies 
in the vicinity of the Fermi level to be so large that the energy positions of 
electrons on hopping sites within some energy strap of width A < A c can 
be arbitrary with respect to the Fermi level. In such a regime, the role of 
quantum mechanical resonant tunneling becomes important leading to the 
universal temperature-independent preexponential factor in the expression 
for the hopping conductivity. 

An attempt to draw a picture of the phononless variable-range hopping 
based on the resonant tunneling has been recently made by Baranovskii and 
Shlimak [9]. They calculated resonant transition probabilities in the regime 
of strong fluctuations of site energies. They also assumed that sites with ener- 
gies in the vicinity of the Fermi level have equilibrium occupations. However, 
the latter assumption seems not consistent with the picture of strong energy 
fluctuations. In the present report we try to derive the variable-range hopping 
transport mechanism which is not based on the equilibrium occupations. 

Our main idea is the following. Due to electron transitions in the "fluc- 
tuators", the energies of sites in the vicinity of the Fermi level, i.e., of sites 
responsible for hopping transport at low temperatures, perform strong energy 
fluctuations within some effective range A. The energy range A increases 
with increasing fluctuation amplitudes. Sites outside this energy range posses 
equilibrium occupations. Transport is due to resonant transitions which oc- 
cur when the fluctuating energy of a filled site within the range A coincides 
with the fluctuating energy of an empty site in this range. The critical hop- 
ping length r c is determined by the solution of the geometrical percolation 
problem on random sites with concentration the latter being the concen- 
tration of sites in the energy strap A in the vicinity of the Fermi level. The 
amplitude of the site energy fluctuations determining A depends on tempera- 
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ture T and it is via A how the temperature influences the conductivity. With 
rising T the width A increases leading to the increase of the concentration 
of effective sites which participate in the resonant-tunneling processes. 
Increasing leads to decreasing r c . This can be well called "the variable- 
range resonant tunneling" [9], because the higher is T, the shorter in space 
are effective hops. Below we give arguments which lead us to such picture of 
hopping transport. 

Let us consider a pair of sites with site 1 initially occupied and site 2 
initially empty. Let the activation energy for electron transition between the 
sites be E and the distance between the sites be R. Transition rate of an 
electron from site 1 to site 2 has the form 

i/i i2 = i/ exp(-2R/a - E/T) (1) 

where z/ is the attempt to escape frequency, a is the localisation length and 
temperature is measured in the units of energy. In the absence of fluctuations 
of site energies the typical transition time t 12 is determined by the condition 
i 12 — I'l l = tc- The fluctuations of the site energies can diminish the 
activation energy of such a pair by some amount SE. Given some spectral 
density of fluctuations, the quantity 5E is time dependent. We assume here 
the spectral density of fluctuations in the form < (SE) 2 >= a 2 (/o//); f < fo 
where / is the frequency and a, fo are some parameters. 

We assume that fluctuations originate mainly due to pair excitations with 
nearly constant density of states [8]. The strongly interacting pairs form a 
sort of "dipolar Coulomb glass". Thus one expects that the 1/ / spectrum 
extends up to very low frequencies, particularly due to the nearly-degenerate 
character of the ground state of the Coulomb glass [10]. The high-frequency 
cut-off fo should be related to the fastest possible modes of the dipolar glass 
related to hops within single pairs independent from the rest of the glass. Cor- 
respondingly, the value of a can be related to the magnitude of fluctuations 
originated due to these single-pair hops. We start from a phenomenological 
picture not specifying the values of a and f and having in mind to return 
to this problem later. 

With an account of such fluctuations the rate equation for electron tran- 
sition within the chosen pair has the form 

dn 2 /dt = u exp {-2R/a — (E + 8E(t))/T) m, (2) 
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where occupation numbers n-i and rt\ have initial values ni(0) = 1; 712(0) = 0. 
This equation has solution 

n 2 (t)= I* dt'v exp{-2R/a- {E + 5E{t'))/T)n 1 {f). (3) 
Jo 

Transition time t is approximately determined by the condition 

f dtV exp (-2R/a - (E + 5E(t'))/T) ~ I. (4) 
Jo 

The integral is obviously dominated by the smallest possible value of 
8E(t). We replace the term SE(t) in the integrand by the mean quadratic 
fluctuation of the activation energy — | < SE(t) > |. Using Eqs.(2), (5), one 
obtains for the typical transition time t 

a(ln(t/ )) 1/2 ~Tln(t c /t). (5) 

Of course, the solution of this equation is related to the typical transition 
time only in the case that the latter is much longer than the time of the pure 
tunnelling without activation th = Vq 1 exp(2i?/ a). Combining this restriction 
with the formal solution of Eq. (6) one obtains 

ln(f/ ) = max(ln(Wo);ln(t c / ) - [(a 2 /2T 2 ) 2 + (a/T) 2 ln(tj ] 1/2 + (a 2 /2T 2 ) 2 ) . 

(6) 

In the case (a 2 /2T 2 ) 2 < (a/T) 2 (ln(t c / )) 1/2 this result reflects a decrease 
of the effective transition time t with respect to its value t c in the absence 
of fluctuations caused by the reduction of the activation energy. In that 
sense our result is similar to that obtained by Ortuno and Pollak [6]. If 
the parameters of a system under study correspond to the opposite case 
a 2 /AT 2 > ln(t c / ) one can use a series expansion in Eq.(6) which leads to 
the expression 

ln(t/ ) = max (\n(t h f ); (\n(t c f )) 2 T 2 /a 2 ) . (7) 

Eq.(7) shows that the transition time in the considered pair corresponds to 
the maximal of the two times: the time of pure tunnelling and the time 
necessary to nullify the activation energy for the transition. This scenario is 
different from the standard VRH where the logarithms of the tunneling term 
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and the activation terms should be compared rather than the corresponding 
times. We assume below that it is condition (7) that governs transition times 
in the pairs essential for hopping transport in the fluctuation regime. 

In order to evaluate the hopping conductivity using given transition rates, 
one has to formulate the binding criterion for construction of the percolation 
cluster. Let us define the binding criterion by the condition that the time of 
the pure tunneling th in a pair is of the order of the time necessary to eliminate 
the activation energy E of the pair. Herewith we refuse to estimate numerical 
constants in the exponential terms for the hopping conductivity, though we 
believe that our criterion gives a correct set of parameters in the exponents 
in analogy to the standard variable-range hopping approach. Our criterion 
corresponds to the optimization procedure that leads to approximate equality 
of two terms in the right-hand side of Eq.(7). 

In order to determine the width A of the energy strip in which electron 
transitions are responsible for transport, let us consider the density of states 
(DOS) function. We assume that DOS at T = is determined by the soft 
Coulomb gap with quadratic energy dependence [11]: 

9(e) = 9oe\ (8) 

where e is the site energy measured with respect to the Fermi level. Its 
integral up to some energy e determines the concentration of cites N £ = 
(2/3)go£ 3 with energies less than e. Thus one estimates as R £ = c/V" 1 / 3 the 
critical distance between the sites within the percolation cluster involving 
the sites from this region of energies (c ~ 0.87, [11]). 

Using this estimate for the first term in the right-hand side of Eq.(7) 
and neglecting logarithmic terms, one obtains that energy ( that provides 
equality between two term in the r.h.s. of Eq.(7) is determined by the relation 
]n(thfo) ~ (4/3)c^ " 1 ^ 3 A _1 /a. In the same way one estimates for the second 
term in the right-hand side of Eq.(7) \n(t c fo) ~ A/T. Thus from the binding 
criterion one has 

A ~ (4c ff o 1/3 « 2 /3a) 1 / 3 . (9) 

It is worth noting that except of a numerical factor the same estimate is valid 
for 2D case. 

According to the transport picture formulated above, the hopping conduc- 
tivity is determined by the solution of the geometrical percolation problem on 
randomly distributed sites with concentration The width of the energy 
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strip A depends on temperature via the temperature dependence of param- 
eter a. We will discuss this dependence later. Before doing so we would like 
to justify our transport picture by showing that in contrary to standard VRH 
approaches, in the suggested transport mechanism the occupation numbers 
of sites which form the transport path are independent on temperature and 
hence these occupation numbers are not essential for the calculation of the 
temperature dependence of the hopping conductivity. 

Taking into account the explicit time dependence of the site energies one 
can write the rate equation for the occupation number of a site i in the form 

df = E (^(^)^(OK(l - m) ' W hJ (e t (t),e,(t))n t (l - nj)) (10) 

where for hopping probabilities one has max W it j = max Wj ti = v exp(—2R it j/a). 
The time dependence of the probabilities does not allow to apply a simple 
detailed balance considerations. Moreover, the 1/f spectrum of the fluctu- 
ations leads to the dependence \5e\ oc ln 1//2 (£_/o) which does not allow any 
time averaging between the hopping events. This is a clear manifestation 
of non-ergodic behavior. The electron hop is prepared by the surrounding 
which implies memory effects instead of ergodicity. 

Energy fluctuations during the typical hopping time do not allow to spec- 
ify any ordering of states with respect to their energies within the energy strip 
A in the vicinity of the Fermi level. As a result, occupation numbers of those 
states do not correspond to the Fermi distribution being rather arbitrary. For 
energies outside the range A , fluctuations cannot suppress the energy or- 
dering of the sites. According to the very idea of the above estimation of the 
hopping time £a, the latter corresponds to the relaxation time of a charge in 
the energy strip of the width A (determined by Eq.(10)) and thus it char- 
acterizes the largest possible time for the charge transfer within the system. 
Having this in mind, one expects that the ergodicity and detailed balance 
considerations are restored at time scales larger than and correspondingly 
for energy scales larger than A. Thus one obtains for the time-averaged 
occupation number of a state i 

J2(< W{ej,£i >< n{£j) > (1- < n(e,i) >)- < Wie^Ej) > ■ 
j 

< n(si) > (1- < n{Ej) >)) = (11) 
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where <> denotes the time average. For the probabilities < W(ei,£j) > one 
obviously has d < W > /d < e^j > \ ei}j ^o << W > /A due to averaging over 
instantaneous values of the energies within the relevant interval ~ A. Thus 
the the averaged occupation numbers < n« > near the Fermi level correspond 
to a smeared Fermi distribution with effective temperature of the order of 
A. For energies |e| >> A the density of states is larger, the tunneling time 
is smaller and the distribution tends to the thermal one. 

The fact that site occupation numbers near the Fermi level are arbi- 
trary justifies the procedure applied above where we neglected the occupa- 
tion number factors restricting ourselves to the diffusion consideration. It is 
also worth noting that according to the above picture, the arguments leading 
to the Coulomb gap [8] do not hold anymore implying that the Coulomb gap 
is suppressed in the energy region of the width A in the vicinity of the Fermi 
level. 

In our model, hopping transitions correspond to processes in which activa- 
tion energy is eliminated by fluctuations of the site energies due to correlated 
electron transitions in the surrounding pairs. In such a case, special atten- 
tion should be devoted to resonant phononless electron transitions. Until 
now, we mainly discussed a preparation of the system of two sites to the 
electron hop rather than the hop itself. One should realize that fluctuations 
with spectral density 1/f contain contributions of all possible modes with 
frequencies in the range between f and 1/t, where t is the hopping time. 
Function ln^ 2 (t/ ) represents a sort of "envelope function" for the fluctua- 
tion amplitudes and activation energy for a hop can vanish many times before 
the hop really occurs. For a pair of hopping sites 1 and 2 with energies E\ 
and e 2 , the routine perturbation theory [13] under adiabatic conditions gives 
the probability p of the resonant transition during a "fluctuation period" 
p ~ P/hKdE/dt)]" 1 , where / is the overlap integral / = J exp(— 2R/a) and 
E = si — e 2 . Making use of the random character of fluctuations and hence 
summing the probabilities for different "periods", one obtains for transition 
rate between sites 1 and 2 

W ~ I 2 /Eh. (12) 

Assume that our sites 1 and 2 form the effective resistor on the percolation 
cluster. Under applied external voltage V, dc current through the chosen 
resistor is equal to 

j = eW < ni(ei(t) + eV, t) - n 2 (e 2 (t), t) >~ e 2 Wd < n(e) > /de, (13) 
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where W is determined by Eq.(13). Taking into account the estimate d < 
n > /de ~ 1/A we come to the following expression for the hopping current 



The distance R in the exponent of the overlap integral / is controlled by the 
geometrical binding criterion discussed above. 

Preexponential term Iq of the overlap integral depends on the character of 
the hopping sites. For sites with the hydrogen-like Coulomb potentials Iq = 
(2e 2 /3na)(r /a), r being characteristic hopping length, k being the dielectric 
constant. In this case Io/A ~ (r/a) 2 >> 1. Different situation can appear for 
sites with screened Coulomb potential. In the case of impurity screening I Q ~ 
e 2 / 'fix ~ A. Using this estimate for I and the estimate for the exponential 
term in the overlap integral from Eqs.(7),(9) one obtains for the hopping 
conductance 



It can look really curious that we pretend to determine preexponential factor 
in the expression for hoping conductivity, while the very derivation of the 
exponent based on the order-of-magnitude binding criterion cannot provide 
us with a correct numerical factor in the exponent. However, our derivation 
of the preexponential factor in Eq.(15) seems to represent just a revision of 
the relation between diffusivity and mobility (Einstein relation) in the non- 
ergodic situation caused by correlated hops with 1/ / spectrum. The crucial 
estimate for the above derivation is the approximate equality between Iq 
and A, both being controlled by the Coulomb energy. The situation is to 
some extent similar to the standard band conduction via extended states in 
a Fermi system with a single characteristic energy equal to the Fermi energy. 

In order to derive the values of T' and (3 it is necessary to specify the mi- 
croscopic picture of the fluctuations which lead to the resonant transitions. 
Following ideas of Ref.2, we assume that the crucial parameter of the fluc- 
tuation spectrum a can be related to the fluctuation potential caused on a 
hopping site by the nearest soft pair with the activation energy ~ T. This 
gives the estimate a = A(AcT 2 )l/3, where A is a numerical coefficient [2, 
12]. Substituting this expression for a into Eq.(15), one obtains 



j ~ (e 2 /h)VI 2 /A 2 . 



(14) 



G ~ (e 2 /h) exp(-(T'/T) /3 ; (T'/Tf = (Ac/Sg^aa) 2 / 3 



(15) 



(3 = 4/9, 



r = (4c/3Aa^ 1/3 A^ /3 ) 3 / 2 . 



(16) 
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Note that the value /3 = 4/9 can hardly be distingwished experimentally 
from the "standard" value (3 = 1/2 for Efros-Shklovskii law. It can be also 
mentioned that for A > 1 the temperature T' appears to be smaller than the 
characteristic temperature Ti for Efros-Shklovskii law; this fact corresponds 
to stimulation of the hopping by fluctuations (additional effect is related to 
an increase of the pre-exponential mentioned above). 

The above estimate of a was based on the parabolic Coulomb gap (Eq.(8)) 
which is believed to exist in 3D systems. For 2D systems, pair excitations 
are known to be less important [11]. Nevertheless, we believe that our results 
can be also applied to many situations where transport takes place in a 2D 
layer surrounded by a 3D system of impurities. Charge hops within the 3D 
dopant system can give rise to the effective energy fluctuations on transport 
sites leading to the transport mechanism suggested in this report. Therefore 
we think that Eqs.(15),(16) might provide a theoretical explanation for recent 
experimental data that show in some cases a universal prefactor e 2 /h in the 
exponential temperature dependence of the 2D hopping conductivity at low 
temperatures [12]. More study is needed, of course, to clarify this question. 
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